function permuted = permute_signal(signal, N, omega, sigma)
%%
% 函数名：permuted = permute_signal(signal, N, omega, sigma)
% 功  能：对信号进行频谱重排
% 参  数：
%       signal：需要频谱重排的原信号
%       N：信号点数
%       oemga：需要取出信号的长度
%       sigma：频谱重排因子
% 返回值：
%       permuted：重排后的信号
%
% 注：这里按照论文，应该还有个随机相位的乘积项 taf，不使用 taf 好像
% 会造成估计相位出现偏差，但还没验证这一点
%%

    permuted = zeros(omega, 1);
    for i = 0:omega-1
        permuted(i+1) = signal(mod(i * sigma, N) + 1);
    end

end